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Abstract. We present a new geometrical method which 
uses weak gravitational lensing effects around clusters to 
constrain the cosmological parameters ft and A . 
On each background galaxy, a cluster induces a conver- 
gence and a shear terms, which depend on the cluster 
projected potential, and on the cosmological parameters 
(ft, A), through the angular distance ratio D^s/Dos- To 
disentangle the effects of these three quantities, we com- 
pare the relative values of the measured ellipticities for 
each triplet of galaxies located at about the same position 
in the lens plane, but having different color redshifts. The 
simultaneous knowledge of the measured ellipticities and 
photometric redshifts enable to build a purely geometri- 
cal estimator (hereafter the G(ft, A)-estimator) indepen- 
dently of the lens potential. 

More precisely G has the simple form of the discriminant 
of a 3-3 matrix built with the triplet values of D^s/Dos 
and observed ellipticities. G is then averaged on many 
triplets of close-galaxies, giving a global function of (ft, 
A) which converges to zero for the true values of the cos- 
mological parameters. 

The linear form of G regarding the measured ellipticity 
of each galaxy implies that the different noises on G de- 
crease as 1/ y/N, where N is the total number of observed 
distorted galaxies. A calculation and comparison of each 
source of statistical noise is performed . 
The possible systematics are analyzed with a multi-screen 
lensing model in order to estimate the effect of perturba- 
tive potentials on galaxy triplets. Improvements are then 
proposed to minimize these systematics and to optimize 
the statistical signal to noise ratio. 

Simulations are performed with realistic geometry and 
convergence for the lensing clusters and a redshift distribu- 
tion for galaxies similar to what is observed. They lead to 
the encouraging result that a significant constraint on (ft, 
A) can be reached : A^q'j in the case ft + A = 1 or ftig 25 
in the case A = (at a la confidence level). In particular 



the curvature of the Universe can be directly constrained 
and the ft = 0.3 and ft = 1 universes can be separated 
with a 2a confidence level. These constraints would be ob- 
tained from the observations of nearly 100 clusters. This 
corresponds to about 20 nights of VLT observations. The 
method is still better adapted to a large program on the 
NGST. Hence, in complement to the supernovae method, 
the triplet method could in principle clear up the issue of 
the existence and value of the cosmological constant 
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1. Introduction 

Determination of the cosmological parameters of the stan- 
dard cosmological models is one of the great challenge for 
the next ten years. Though these are the main objectives 
of the MAP and Planck Surveyor satellites, considerable 
efforts are devoted to the measurements of the (ft, A) 1 
parameters prior to the launch of these surveyors. In this 
respect, the supernovae search (see the Supernova Cos- 
mological Project, Perlmutter et al. 1998) or gravitational 
lensing surveys (for a review see Mellier 1998) offer the 
best perspectives. In particular, if they are used jointly, 
a reasonable precision can be reached, provided that the 
degeneracy regarding to the determination of parameters 
(ft, A) of each method are orthogonal. This is for instance 
the case for the supernovae experiments and the weak lens- 
ing analysis as presented hereinafter or as produced by 
large scale structures (Van Waerbeke et al 1998). 
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1 is the matter density to critical density ratio. A is the 
ratio between the matter density associated to the cosmological 
constant A and the critical density. A is equivalently defined 
as the ratio A/3Hq, where H is the Hubble constant. 
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The limitations of any of these methods (even Planck mea- 
surements) are the understanding of the systematic bias 
and eventually the control of the large number of free pa- 
rameters attached to each of them. Due to the difficulties 
to handle these issues it is important to diversify the meth- 
ods to measure (fl, A) and to find out some new observa- 
tional tests. In this regard, any new method that controls 
properly its own systematics and that decreases the num- 
ber of sensitive parameters needs careful attention. 
It is well known that gravitational lensing can provide 
purely geometrical tests of the curvature of the Universe. 
Applications of this property to lensing clusters have been 
proposed by Breimer & Sanders (1992) and Link & Pierce 
(1998). They suggest to use giant arcs having different 
redshifts to probe directly the curvature of the Universe. 
Evidently this method can provide the cosmological pa- 
rameters in the simplest way, provided that the modeling 
of the lens is perfectly constrained. Besides, it requires the 
spectroscopic redshifts of at least two different arcs in the 
same clusters, which is not an easy task. As such, this is 
a method which applies to very few clusters. 
Fort, Mellier & Dantel-Fort (1996) focused on a statistical 
approach which explores the magnification bias coupled 
with the redshift distribution of the sources. Fort et al. use 
the shape and the extension of the depletion curves pro- 
duced by the magnification of the galaxies to constrain the 
cosmological parameters and the redshift of the sources si- 
multaneously. The intrinsic degeneracy can be somewhat 
broken if the redshift of a giant arc is known and if the 
number density of high-redshift background galaxies is sig- 
nificant. However, in practice, reliable results need the in- 
vestigation of a significant number of lensing arc-clusters 
in order to improve the statistics, to minimize the system- 
atics (like multiple lens planes) and to explore the sensi- 
tivity to the lens modeling. 

The key issue on the Fort et al. approach is the coupling 
between the cosmological parameters, the redshift of the 
sources and the lens modeling. An attempt to disentangle 
these three quantities has been proposed by Lombardi & 
Bertin (1998), using a method which applies in rich clus- 
ters of galaxies, inside the region where the weak-lensing 
regime is valid. They use the knowledge of the photomet- 
ric redshifts for a joint iterative reconstruction of the mass 
of the cluster and the cosmology. However their iteration 
method assumes that the mass of the deflecting cluster is 
known (or equivalently they assume that the mass-sheet 
degeneracy inherent to the mass reconstruction is broken). 
This assumption is the key of the problem, because it 
means that a perfect correction of the systematic bias on 
the (fi,A) determination due to the systematic errors in 
the mass reconstruction (see section 3.1) can be achieved, 
which is still not presently the case (see Mellier 1998 for a 
comprehensive review). Indeed, for such a curvature test 
one should try to escape to the uncertainties of the poten- 
tial modeling (mass reconstruction)of the cluster. 
The triplet method proposed in this paper can solve this 



problem because it is based on the construction of an 
(il,A)-cstimator independent of the lens potential. Basi- 
cally, it consists in comparing the elliptical shear of 3 
nearby galaxies of different color redshifts probing the 
same part of the potential across the cluster. When the 
galaxies are close enough, their observed ellipticity only 
depends on three unknown parameters, the local conver- 
gence and shear (related to the second order derivatives of 
the projected potential of the cluster) and namely the cos- 
mological parameter (fi, A), through the angular distance 
ratio Dls/ Dos- The use of a triplet of nearby galaxies 
enables to break this degeneracy and provides a local ge- 
ometrical operator only dependent on (f2, A). We chose a 
local linear operator regarding to the observed elliptici- 
ties, and we average it on all set of close-triplets detected 
in many lens planes. Follows a global geometrical estima- 
tor, biased by a noise (coming mainly from the intrinsic 
source ellipticities) which decreases as the inverse of the 
square root of the number of triplets. It means that with 
a large number of lensing clusters one can estimate (fi, A) 
with a reasonable accuracy. This technique and its effi- 
ciency are discussed in the following sections. 
Section 2 reminds rapidly the lensing equations and the 
basic lensing quantities relevant for the paper. Section 3 
shows how to build the geometrical estimator G which 
uses triplets of distorted galaxies. The principle and the 
detailed analysis of our method are also discussed. The 
signal to noise ratio of the method is then derived as well 
as the probability distribution of the cosmological parame- 
ters (fi,A). Although the main objective of this first paper 
is mostly to present the principle of the method, Section 4 
gives an evaluation of the amplitude of several systematic 
biases coming from possible perturbing lenses distributed 
along the line of sight of triplets (galaxies or larger struc- 
tures). Preliminary solutions to these systematic biases 
and ideas of optimizing the method (Selection in the ge- 
ometry of clusters and choice of redshifts) are developed 
in section 5. The method is tested on simulations in Sec- 
tion 6. Finally, we discuss the results and suggest some 
observational strategies in Section 7. 



2. The weak lensing equations 

The lensing properties are determined by the dimension- 
less convergence (the strength of the lens) n and shear (the 
distortion induced by the lens) 7 2 , which both depend on 
the second order derivatives of the two-dimensional pro- 
jected deflecting potential. The lensing effect of a cluster 
on background galaxies can be expressed as an amplifica- 
tion matrix defined in each angular position around the 

2 In the following mathematical notations with bold letters 
refer to complex numbers while usual letters are used for scalars 
or for the norm of the associated complex numbers. The up- 
per * index behind a complex number indicates its conjugated 
element. 
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cluster as (Schneider et al. 1992) 

A -i = A-k-71 -72 

\-j2 1 - K + 7i 

where 

7 = 7i + *72- 



(1) 



(2) 



k is a dimensionless form of the cluster surface mass den- 
sity £ : 



k — Y*/Yj cr it with T, cr n 



Dos 



AttGDql D ls 



(3) 



where Dos, Dol and Dls are the angular diameter dis- 
tances from the observer to the source, from the observer 
to the lens and from the lens to the source respectively. 
For the weak lensing regime the gravitational distortion 
produced by a lensing cluster can be modeled by a trans- 
formation in the ellipticity of the galaxies from the source 
plan (es) to the image plane (e) (see Appendix A) : 



e = (1 - g 2 )e s + g = e s + g with g 



7 



1 - K 



(4) 



e is the complex observed ellipticity, g is the complex re- 
duced shear and e$ is what we will call the complex cor- 
rected source ellipticity. Either in the source or the image 
plan, the ellipticity parameter is defined by : 



e = e e 



2,0 



with e 



1 



1 + r' 



(5) 



r is the axis ratio of the image isophotes and 6 is the 
orientation of the main axis. 

The convergence and the shear both depend on the source 
rcdshift through an absolute lensing factor u a appearing 
in equation (3) : 



w a (z) = 



Dls 
Dos 



(6) 



Adopting the notation of Seitz & Schneider (1997) who 
relate the lensing parameters to the value they would have 
at infinite rcdshift, k and 7 now write : 



K 

1 



w ( z )7o. 
w (z) 



where 



(7) 

(8) 

(9) 



w a (oo) 

Hereafter, u>(z) will be named the lensing factor. This 
is the term which contains the cosmological dependency. 

3. The method 

The behavior of oj(z) with fi in the case Vt + A = 1 (flat 
geometry) and in the case A = 0, is given on figure (1). 
All curves range from (for a source redshift equal to the 



cluster redshift) to 1 (for an infinite source redshift). Their 
main difference is a small change in their convexity. That 
is why the use of numerous triplets of sources at different 
redshifts may disentangle these curves, i.e. provide a con- 
straint on the cosmological parameters. 
The core of the method is to proceed in such a way that 



cj(z) Case f)+A=1 n = 0.2,.. 1 .0 
1 .0 



w(z) Cose A=0 [5 = 0.2,.. 1 . 





2 3 4 



Fig. 1. u) versus redshift in the cases £1 + A = 1 and A = 0, for 
a cluster lying at a redshift 0.4. Different curves correspond 
tofl = 0.2 (bottom curve), 0.4,0.6,0.8 and 1.0 (top curve). 
All curves start from (at redshift 0.4) and converge to 1.0 
when the redshift becomes infinite. The triplet method applies 
to their different relative convexity. 



this constraint is independent on the potential of the clus- 
ter. We have constructed an operator which depends the- 
oretically on (ft, A) and can be computed simply from 
the observed ellipticities of background galaxies as well 
as their photometric redshift. Its main property is to be 
equal to zero when the cosmological parameters are equal 
to the actual ones. We proceed in two steps: first, we build 
such an operator (Gijk) from triplets of close sources, and 
second we average it on many triplets of sources to obtain 
the final geometrical operator G. 

3.1. Construction of Gijk 

For a background galaxy at redshift z$, equation (4) 
rewrites 



e 4 = e SjJ + g° with 

1 - W?Koo 



(10) 

(11) 



where the lower index (i) refers to the redshift Zi and the 
upper index (°) refers to the actual values of the cosmo- 
logical parameters : o ,A o . The second term in equation 
(10), g° represents the part of the image ellipticity that 
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depends on the cosmology. 

Let us now consider a triplet of background neighboring 
galaxies in the image plane and lying at redshifts z%, Zj 
and Zk- The number density of triplets depends on the 
deepness of the observations. For instance, up to R = 26 
mag., we expect a mean density of 50 sources arcmin~ 2 , 
that is about 4 sources inside a circle of radius 10 arcsec- 
onds. In the following, we thus consider that each galaxy 
of the triplet is distorted by the same local potential i.e. 

and 7^ are the same for the three galaxies. The bias 
induced by this approximation will be discussed in sec- 
tion 3.4. The triplet of galaxies gives a triplet of equations 
(11) respectively indexed by i, j and k from which we can 
derive a final equation independent both on /-too and 7^ . 
This equation writes simply as the zero of a 3-3 determi- 
nant : 



1 < 

1 w? uffig? 
1 < u° k 9?9° 



(12) 



The first term of this equation can now be formally gen- 
eralized to a complex operator Gijk of (f2,A) built from 
the complex measured ellipticities of the three galaxies : 



G ijk (fl,A) 



1 uji bJie^e 



3 c k 



1 uj uJje k e* 
1 uj k ^fc£ 4 e* 



(13) 



The dependency in (fl, A) is contained in each term 
u> s (Ct, A) (s = i, j or k) defined by equation (9) for z = z s . 
djk is more explicitly the sum of two (ft, A) functions: 



G 



ijk 



which is equal to zero for the actual values of the 



cosmological parameters, and a complex noise A^: 



ijk — ^ijk 



G 
where 



+ N 



/-imam 
^ ijk 



1 OJi 0Ji/0J° 
1 OJj OJj/oj] 
1 LU k 0J k /bjl 



7^. 



(14) 



(15) 



At this point it is important to stress that the Koo contri- 
bution cannot be neglected (see equation (11)). Indeed, if 
we consider for instance a 0.1 variation of the cosmological 
parameters (along the gradient of w(fi,A)), the resulting 
relative variation of the term is about 1% which is 
about ten times smaller than the relative variation due to 
the 1 — ujk term, in equation (11). That is why both con- 
tributions from the local convergence («oo) and the local 
shear (700) of the cluster potential must be taken in ac- 
count. 

Lombardi & Bcrtin (1998) proposed to reconstruct jointly 
the shear, the convergence and (CI, A) in the weak lens- 
ing area, with an iteration method based on the equation 
(11). Their method seems to converge relatively rapidly 
with a small number of clusters but it seems that for their 



simulations they implicitly assume that the mass of the 
cluster is known. However, one can see that equation (11) 
is invariant when replacing 700 by 0:700 and 1 — wk m by 
1 — aujKoc (a is a constant). This expresses the mass- 
sheet degeneracy problem which implies that the total 
mass of the lensing cluster is uncertain. Indeed, despite 
the numerous suggestions which have been proposed in 
order to solve this issue (Seitz & Schneider, 1997), for 
the moment mass reconstruction techniques cannot dis- 
regard systematic errors analysis (see Mellier 1998 for a 
comprehensive review) . As an example, a 20% systematic 
bias on the determination of the total mass of the cluster 
(or equivalently a 20% systematic on the mean value of 
Koo) is equivalent to a systematic bias larger than 0.2 on 
the value of the cosmological parameters (when compared 
to the 1% contribution from the lensing factor mentioned 
above). Therefore, the knowledge of the lens potential is 
a critical strong assumption. 

In the triplet method it is possible to constrain the cos- 
mological parameters regardless the potential of the lens. 
No assumption is made in order to relate the values of the 
local shear and the local convergence. Hence they are con- 
sidered as independent parameters. In order to construct 
a G operator which only depends on (f2, A), we then need 
three local equations relating Kqq, 700 and (f2,A) to can- 
cel the potential dependency. This is achieved with the 
measured ellipticity equations (10) applied to triplets of 
close galaxies at different redshifts. Besides, the form of 
the G operator becomes unique and must have the formal 
expression given in equation (13) if we want it linear with 
respect to the ellipticities provided by the observations. 
To sum up, the use of triplets of galaxies through the op- 
erator G is the simplest way to build a pure geometrical 
operator which drops both the Koo and 700 dependencies 
and keeps linear regarding to the ellipticities. 

The statistical noise TV will be more explicitly calcu- 
lated in section 3.4. For the moment it is just important to 
understand that the probability distribution of this noise 
(real and imaginary parts) regarding the different triplets 
of galaxies is a random law centered on since the lin- 
ear construction of G makes the different sources of noise 
(mainly the intrinsic source ellipticity) be randomly dis- 
tributed around 0. 

The above formula does not take in account the systemat- 
ics (an effect of galaxy-galaxy lensing, a presence of back- 
ground structures) that will be studied further in section 
4. 



3.2. Construction of G 

Before averaging Gjjfc(Q,A) on many triplets (that is to 
say on all the triplets done with galaxies contained within 
a given radius), let us stress that it makes sense. Indeed, by 
averaging (we only consider the real part of the complex 
noise) on many triplets the noise contribution decreases as 
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1 / VN, with N the number of background galaxies (and 
not the number of triplets, since many triplets are redun- 
dant). While the noise is vanishing, the method can pro- 
vide valuable constraints as far as cosmological operator 
G™% m does not vanish too, or in other words if its behav- 
ior with SI and A is the same (not random) whatever the 
triplet is. Fortunately this is the case provided that the 
three redshifts of all the triplets are ordered similarly : for 
example Zi < Zj < z k . 

Consequently we can derive G from an average of Gyfc on 
all the ordered triplets ijk : 



G(Sl,A) =< djk >{i<j<k) 



(16) 



According to equation (14), G is also composed of two 
terms : the first one G mam which has the same properties 
as G^f™, and a Gaussian noise GN which decreases as 



(17) 



l/VN, 

R e (G)(Sl,A) = G mam (Sl, A) + Re(GN) 
where 



G" 



and 



1 LUi UJi/uJ° 
1 Wj Uj/uj] 
1 UJ k uj k /uj° 



(ijk) 



Re(GN) OC l/VN , 



(18) 



(19) 



where R e denotes the real part of the complex quantities. 
< 7^ > is the square of 700 averaged on the weak lensing 
area. We consider the real part of G (Re(G)) because as 
N, the noise GN is complex. 

By construction G mam is equal to zero at the position (Sl a , 
A ). We thus can write the following equation (biased by 
the presence of the Gaussian noise) which sums up the 
method : (Sl , A a ) is a solution of 



G(Sl,A) = 0. 



(20) 



In order to see if the triplet method can be effective from 
the observational point of view two questions have to be 
addressed: 

1. how is the G-operator degenerated in (SI, A)? 

2. how many clusters are necessary to cancel the noise 
contribution with respect to the cosmological informa- 
tion contained in G mam ? 

3.3. The main cosmological term: G main (Sl, A) 

Figure (2) gives the contours of G mam / < 7^ > in the 
(SI, A) plan, for the redshift distribution (42), as taken in 
the simulations (see section 6). From this graph we see 
that the method is degenerated in (SI, A). The degeneracy 
is parallel to the G-contours. To break this degeneracy one 
has either to make a theoretical assumption (for example 



SI + A = 1 ) or to add another experimental degener- 
ated constraint which contours are as orthogonal to the 
G-contours as possible, like for example high-redshift su- 
pernovae (Perlmutter et al. 1998). Considering the mean 
orientation of its degeneracy, the triplet method can also 
be seen as a direct measure of the curvature of the Uni- 
verse 1 — SI — A. 

From figure (2) we also get quantitatively the variations of 
G. It shows that, for an accuracy of about 10% on the cos- 
mological parameters (along the gradient of G), we must 
obtain a precision of about 10~ 4 < 7^ > on G. This rate 
of variation has to be compared with the noise on G. 




Fig. 2. Contours of G main (SI, A) / < 7^ > in the cases 
(SI, A) = (0.3, 0.7) (top panel) and (0.3, 0) (bottom panel). The 
right part of the figure gives the restriction of G mam to the 
particular domains Q + A — 1 (top) and A = (bottom). 



3.4. Noise on G 

The complex noise is produced by four sources: 

1. noise from the corrected source ellipticities es,ij,k, 

2. the errors propagation on the measured ellipticities 
^ e i,j,k (it behaves similarly as the previous ones), 

3. the three sources do not have the same 700 and we 
thus have to consider A'-foo^j^ (in other words , each 
source, though close to each other, do not exactly cross 
the potential at the same position), 

4. the photometric redshifts are not the true redshifts and 
lead to shifts Auij^k on the lensing factors. 

From equations (13) it is clear that, due to the linearity 
of the 3-3 determinant and of the averaging on all the 
triplets, the final noise is linear regarding to each indi- 
vidual term : it is composed to first order of linear com- 
bination of terms (like es,i), and to second order of lin- 
ear combination of crossed-terms (like es^.Awj). These 
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crossed-terms cannot introduce any systematic bias on G 
and thus can be neglected. 

The following equations give the four different noise con- 
tributions: GN es A , due to the corrected source ellipticity, 
GN Aei , due to the error on the ellipticity, GN AuJi due to 
the approximate redshift and GN Aloo ( due to the poten- 
tial, 



GN Aei 
GN Auli 



-<Too > 



-<loo > 



-<Too > 



-<Too > 



< - uj k )e* S i /j* (21) 

+u?u k {wj - w l )e s , i /7 > 

<w J -wg(o; i -w fc )Ae,V7* (22) 

< Au>i(u)jUJ^(uJh - Uj) (23) 

<w j wg(wi-Wfc)A7j/7* (24) 



The terms of index j and k can be easily derived from these 
equations by cyclic permutations of the (ijk) triplet. 
It is worth noting that even if the bias, as G mam , is 
a function of the cosmological parameters, we can ne- 
glect this dependency. Indeed, for a given shift of (ft, A) 
around (ft ,A ) the variations of G mam and GN verify 
AR e (GN)/AG main cx 1/VN. Besides, provided that N 
is large enough, the probability distribution of the noise 
R e (GN) can be considered as a Gaussian law centered on 
(as the sum of a large number of nearly Gaussian laws). 
So far, we have shown that the variance of this noise de- 
creases as 1/VN because of the redundancy of triplets. 
This applies for G7Vj S i , GN Aei and GN AuJi but not for 
GN Alooi . Indeed, for a galaxy i included in two different 
triplets the associated term A7^ is different in each triplet. 
That is why the variance of the fourth noise GN Ajoo i van- 
ishes in 1/ y/3N tr , where N tr is the total number of triplets 
: Ntr ~ 3iV (see the remark of section 3.1). 
The variances of the four noises behave approximately as 
follows (under the conditions that will be detailed in sec- 
tion 6) : 



GN- es . 



GN Au 
GN Aulz 
GN Al<Xti 



0.08 



0.02 



0.04 



0.02 



N 



N 

<llo> 



N 



N 



(25) 
(26) 
(27) 
(28) 



3.5. Resulting signal to noise ratio 



Let us establish the relation between the signal to noise ra- 
tio of the method and the number of background galaxies 
(or the number of considered clusters). In the following we 
will take only 1000 galaxies per cluster (see section 6). If 



we define the signal as a variation of G mam along its gra- 
dient and a on as the variance of the statistical noise (see 
section 3.4), the signal to noise ratio can thus be written 
as : 



S/N 



AG" 



A. 



Q,A 



main 



0~ GN 



0~GN 



(29) 



The ||...|| notation indicates the norm of a vector. Ao.a is 
the accuracy on the cosmological parameters. Figure (3) 
plots the variation of Ao.a as a function of the number of 
clusters (with an observed number density of background 
galaxies equal to 50 by square-arcminute) for a lc confi- 
dence level. This plot shows that an accuracy of 0.1 can 
be reached (at a la level) with about 1000 clusters, or 
similarly an accuracy of about 0.3 can be reached (at a la 
level) with 100 clusters. 

It shows that with a good seeing one could in principle test 
the existence of a cosmological constant from VLT obser- 
vations with about 100 clusters. Far better it should be 
possible to measure the curvature of the Universe with a 
reasonable accuracy from NGST observations (the num- 
ber density of background galaxies is multiplied by 10). 
These prospects of course require that the systematics are 
previously corrected. These results can be refined by the 



Accuracy on (fi,A) 



400 600 800 1000 

Number of clusters 



Fig. 3. Accuracy on the cosmological parameters versus the 
number of clusters (with 1000 galaxies per cluster) at a la con- 
fidence level. The plain line extrapolates the results obtained 
from simulations with an inverse square root function. Error 
bars account for the number of simulations. We can read that 
a 0.3 accuracy (at la) can be achieved with the observations 
of 100 clusters or a 0.1 accuracy (at la) can be achieved with 
the observations of 1000 clusters. Equivalently a 0.2 accuracy 
(at 2a) can be achieved with the observations of 1000 clusters. 



introduction of a probability distribution. 
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3.6. The (fi,A) probability distribution 

Thanks to the Gaussian behavior of the noise we can de- 
rive from the operator G(fl, A) a probability distribution 
A) for the cosmological parameters : 

pn , A(aA)sexp (_l5W)!f), ,30) 

which has just to be normalized on the considered domain 
of (Q, A), (for example Q + A = 1). 

The operator G(fl, A) is directly derived from the obser- 
vational data. On the other hand, <jgn is obtained from 
simulations applied on the data in the following way. Once 
the function G(£1,A) is obtained from the observation of 
(e, z), it is possible to add noise to the data by adding 
random ellipticities (es) and Az random errors to the red- 
shifts z : the variance of the noise GN induced on G is 
then directly ctgn- 

4. Analysis of systematics 

In this section we discuss the systematics which could po- 
tentially reduce significantly the efficiency of the method 
from an observational point of view. Though we give es- 
timations of their amplitude, it is worth noting that a re- 
liable quantitative estimate of their impact on the triplet 
method will demand additional work (mainly simulations 
and calculations taking in account the perturbation com- 
ing from large-scale structures. This effect will be studied 
in paper II). 

Two main systematics have been identified: 

1. a systematic bias on G produced by to the non- 
symmetry of the term Aoj. It can be easily corrected 
(see 4.1). 

2. a contamination produced by the existence of back- 
ground structures which play the role of perturbing 
lenses. They can be the background galaxies them- 
selves (which induce galaxy-galaxy lensing) and pos- 
sibly any condensation of mass (clusters or large-scale 
structures) located on the line of sight of the lens- 
ing cluster. In the following sections (section 4.2 and 
4.3), we focus on two extreme cases: the systematic 
produced by the galaxies of each triplet themselves, 
and those generated by large-scale structures. This 
study will be performed using multi-lensing models 
(see Kovncr 1987). Appendix B gives the calculation 
of the measured ellipticities in a multi-lensing model. 

4-1. Non symmetry of the term Auj 

The method supposes that we know the photometric red- 
shifts of background galaxies as well as the probability 
distribution p& z {Az) of the error between the photomet- 
ric rcdshift and the true one. Even if p Az {Az) is symmet- 
ric and centered on 0, which is not necessarily the case, 



the probability distribution of the resulting shift on uj, 
PAuj(Auj) may be non symmetric and the mean value of 
Au) may be different from zero, introducing a systematic 
in equation (23). Fortunately, this effect is easy to correct. 
Since the redshift distribution n(z) is known the system- 
atic can be exactly balanced by replacing in the definition 
of Gijk u) by uj — Auj, where 

Aoj(z) = J uj(z + Az)n(z + Az)p Az (Az)dAz . (31) 

4-2. Potential perturbation due to galaxy-galaxy lensing 

Since we consider triplets of galaxies having different red- 
shifts but which are very close together in the image plane, 
the most distant can potentially be distorted by the lens- 
ing effect of the others. In particular, in each triplet k) 
(ordered with increasing redshifts) the galaxy i is a per- 
turbation of the potential for the sources j and k as well 
as the galaxy j for the source k. It is worth noting that 
i is not the only perturbing lens for the source j, and i 
and j are not the only one for the source k. For instance, 
each galaxy may be itself part of a group of galaxies which 
produces its own additional lensing effect. Nevertheless, in 
order to get a rough estimate of the amplitude of this kind 
of systematics, we will assume to first approximation that 
the galaxy-galaxy lensing produced within each triplet, is 
the dominant contribution. 

The galaxy-galaxy lensing increases as the angular dis- 
tance between the galaxy-lens and the galaxy-source de- 
creases. It also depends on their relative redshifts: that is 
why it induces a systematic bias on the construction of 
G. In Appendix B we compute the amplitude of the cor- 
rections on the measured ellipticities once several lenses 
along the line of sight are taken into account. From Equa- 
tions (62), (63) and (66) of Appendix B, one can see that 
these perturbations produce two kinds of taxes: 

1. A purely additive term, g p . It can be seen as the lin- 
ear contribution from the shear, as if the shears of each 
lenses were simply added to each other without cou- 
pling considerations. 

2. A a multiplicative term of the form (l — (1 — c) k p ) 
which modifies the main cosmological term g. It ac- 
counts for couplings between the main and additional 
lenses. 

4.2.1. The additive linear term. 

Since the three galaxies of each triplet are randomly dis- 
tributed and uncorrelated, this term induces a noise rather 
than a systematic bias. It decreases as l/y/3N, like the 
noise GN Ajoo for instance. Besides, the shear induced 
by a galaxy-lcnsed can be controlled (if we cancel the 
triplets where the angular distance between two of the 
three sources is too small: typically lower then a few arc- 
second for typical Einstein radius of galaxies. This means 
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an induced shear inferior to about 2%). Therefore the 
amplitude of this additive term can be bounded to less 
than 10% of the effect coming from the principal statisti- 
cal noise GN es i5 and thus can be neglected. 



4.2.2. The multiplicative coupling term. 

The effect of coupling between two lenses depends on their 
relative redshifts through angular distances ratio, hence 
this term really induces a systematic bias on G. To calcu- 
late it we have to think about the meaning of the addi- 
tional convergence, n p ', used in the multi-lens screen ap- 
proach of the annex B. This convergence associated to the 
source j and produced by the galaxy-lens i will be noted 
k 1 ' ; the convergences associated to the source k and pro- 
duced by the galaxy-lenses i and j will be noted n l,k and 
K j,k_ -yy e cannot directly compare their value since they 
depend on the redshifts Zi,Zj,Zk- So, in order to model 
this systematic bias we associate to each of these conver- 
gences an absolute one k gal independent of the redshifts 
Zi, Zj and Zk (obviously k gal should depend on the angular 
distance between the galaxy-source and the galaxy lens. 
We consider here a mean value averaged over the area 
around each galaxy-lens where we search for a source. For 
the definition of this area, see section 5.2): 



Dol Doj 

Dot Dik ~ ga l 
— f% 

Dol Dok 
Doj D 3 k ~ ga i 



Dol Dok 



(32) 
(33) 
(34) 



The Aq ^ L shift (along the gradient of G) on the cosmo- 
logical parameters determination and due to the galaxy- 
galaxy lensing perturbation effects is then : 



Ldi 

1 T& 





Em (Dm _ Dm i 

Dol \ D 0j D Lj > 
Di 



I Dp % _ Phi \ I Djk I Doj _ DLj \ 
V D ok D,, h > T Dol 1 D„ k D,, k > 



VO ,A 



1 LJi UJi/uj° 
1 Uj UJj/uj" 
1 UJ k UJk/^l 



"(35) 



where the brackets (...) denote the average over all the 
ordered triplets. 

Using the redshift distribution (42) and considering the 
mean value of the convergence created by a galaxy as 
about 2% (see section 5.2), one gets : Ag G A L = 10.2k gal rj 
0.2. 

This is only an indicative value. However, with 
galaxy/galaxy lensing measurements, we are in principle 
able to account for the effect of the potential of each back- 
ground galaxy. This systematic can be calculated and cor- 
rected. Therefore, if the knowledge of the mean potential 



of galaxies can be reached with about a w 20% accuracy 
(we could use Faber-Jackson models or wait for the incom- 
ing results of galaxy-galaxy lensing studies (see Schnei- 
der et al. 1996)), the remaining systematic would then be 
about Ag G A L w OM/VNclust. 

This correction does not take in account the possibility 
that a fraction of the galaxies are embedded in groups 
which could enhance the contamination. There arc two 
possibilities to solve this issue. Since we know the photo- 
metric redshifts of the galaxies we can discriminate com- 
pact groups of galaxies and remove the corresponding 
triplets from the sample. Or, if the galaxy-galaxy lensing 
results are able to account for these effects with a reason- 
able accuracy, then we can as previously calculate A^^ L . 
In conclusion this effect is non-negligible but can be accu- 
rately estimated and corrected, thanks to the forthcoming 
observations of galaxy-galaxy lensing on blank fields. 



4-3. Potential perturbation due to background lensing 
structures 

Let us consider the perturbation due to background struc- 
tures. Before giving an estimate of it for the case of matter 
structures integrated along the line of sight, we calculate 
the perturbation due to a single lens plane (containing an 
over dense region like a cluster or even an under dense 
region). Since this second lens may be located anywhere 
in redshift then it differently affects the measured elliptic- 
ities of the sources k) in each triplet, it can induce a 
systematic on the value of G. Once again and still using 
the calculations done in annex B (see equation (66)), this 
systematic can be split in two parts due to an additive 
term (corresponding to the linear contribution of the per- 
turbative lens) and a multiplicative term (corresponding 
to the coupling between the main and the perturbative 
lenses) . 



4.3.1. The additive linear term. 

Like in the above section, we associate an absolute conver- 
gence and shear 5k and 5j to the additional convergence 
and shear Sk P ' S and <57 P ' S coming from the lensing effect 
of the structure on a source S (S can be i, j or k): 



5k p > s = 



Dql p Dl p s 
Dol Dos 

Dql f Dl p s 
Dol Dqs 



8k 



(36) 
(37) 



where the superscript p denotes the perturbative term. 
We assume that the convergence and the shear are con- 
stant all over the image. The A^ B A S shift (along the gra- 
dient of G) on the cosmological parameters determina- 
tion and due to the linear lensing effect of the background 
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structure is then: 

1 Wi 
1 Uj 
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For a condensation of mass located at redshift 1 and 
with the redshift distribution (43) this systematic is : 

A^ s = 1.8 (jf^j ■ To generalize this systematics to 

large scale structures, the above value of A^^ s has to be 
integrated along the line of sight. Here we do not perform 
exactly this calculation (which needs the introduction of 
cosmological scenario and the non linear evolution of per- 
turbations) but only give an estimation of it using the 
already known results from the studies of lensing by large 
scale structures (see the review by Mellier 1998). Calcula- 
tions from the non linear evolutions of the power spectrum 
in the one minute angular scale predict a polarization of 
about 3%. From this value we can make a rough estimate 
of the shift on (Cl, A) (using the same conditions as will 
be used in the simulations, section 6) : A^ s w 0.03. 
This is only an indicative value, however this systematic 
can be estimated more accurately using the simulations of 
the non-linear evolution of perturbations. 

4.3.2. The multiplicative coupling term. 

The A^' s shift (along the gradient of G) on the cosmo- 
logical parameters determination and due to the coupling 
effect of a background perturbing structure (with the main 
lens) is calculated as in 4.2.2 but is simpler since there is 
only one perturbating lens to which we associate an ab- 
solute convergence and shear (see equation (66)). It leads 
to: 
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(39) 



The same discussion as above applies. For a condensa- 
tion of mass located at redshift 1 and with the redshift 
distribution (43) this systematic is : A^ s = 16.75k. To 
generalize this systematic to large scale structures, the 
above value of Aq ^ s has to be integrated along the line 
of sight. We use again a polarization of about 3%. Here 
the difference is that this systematic effect behaves like 
a statistic noise because the SR can be either positive or 
negative. Hence, when the triplet method is applied on 
many different clusters, we expect the averaged value of 



this systematic to behave roughly as the inverse of the 
square root of the number of considered clusters. As an 
example, for 100 clusters, we expect this systematic to be 
about Ag^ s w 0.04. 

Contrary to the other two systematics (Aq^ l and Aq^ s ), 
this one can not be corrected because of its random be- 
havior. However its effect on the (CI, A) determination de- 
creases with the number of considered clusters. 

4.3.3. Effect of a background cluster 

In the last two sections we discussed the non-linear evo- 
lution of background structures in a scale of about one 
arcminute to estimate the resulting systematics on the 
triplet method. It may be also important to take in ac- 
count (statistically) the effect of background clusters. In 
particular, the accidental presence of another cluster lying 
along the line of sight of the main lensing-cluster could be 
a serious artefact. In the case of usual clusters of galaxies, 
these biased targets can be easily removed from the sample 
by using photometric redshift informations. However, if it 
happens that dark clusters do really exist, then their pres- 
ence could be extremely difficult to detect. Their impact 
on the statistics depend on the number density of such 
systems, if any. Since we do not have any evidence from 
shear map that dark clusters exist, this is a somewhat an 
academic problem. However, such structures could be re- 
vealed by using a method like the aperture mass density. 
In principle, this technique is able to detect a dark halo, 
provided that their velocity dispersion is larger than 400 
km/s (see Schneider et al. 1995). 

A detailed quantitative estimate of these effects require 
simulations and also additional observations in order to 
constrain the fraction of lensing clusters which is contam- 
inated by another cluster along the line of sight. This in- 
vestigation is beyond the scope of this paper, but should 
be addressed in a forthcoming paper. 

4-4- Outcome 

Estimations of the various systematic biases are given in 
the following table, for a principal lens at redshift 0.4, and 
a redshift distribution similar to (42): 
Next section will give qualitative solutions to deal with 
part of these systematics and to increase the signal to noise 
ratio of the method. 



5. Optimization of the method. 

So far, the method has been presented in a structure as 
general as possible, with no particular choice for the ge- 
ometry of clusters, neither for the redshifts nor for the 
distances between galaxies inside triplets. This section dis- 
cusses qualitatively these degrees of freedom in order to 
increase the signal and decrease the bias in the same 

time. 
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A GGL 


~ o.2 


rorrprtpd to pa 04/a//V > 


\LBS 


w 0.03 


corrected 


A CBS 


« 0.4/ ViV c i ust 


not corrected 



Table 1. The systematic biases on the determination of (fi, A). 
Aq ^ l is due to galaxy-galaxy lensing effects. Aq B a is due 
to the linear effect of the distortion induced by background 
structures. Aq^ S comes from the coupling of these background 
structures with the main lens. The corrections of Aq G a l and 
Aq5l S require to adjust the method with a modeling of the 
potential of galaxies and large scale structures. Such modeling 
compared with ray tracing simulations will be performed in 
paper II. Aq^ s can not be corrected because of its random 
behavior. 



5.1. Choosing clusters with symmetrical geometry 

The method as explained so far can be applied to every 
cluster regardless its geometry. However new cluster sam- 
ples, like the one obtained by the X-ray satellite XMM, 
will be useful to select only those with a symmetrical ge- 
ometry, since any geometrical information can be used 
to decrease the systematic bias noticed due to the linear 
perturbation effect of background structures (see section 
4.3.1). This bias is due to scalar products as g p ' l g p 'i* ap- 
pearing from the construction of Gijk in the determinant: 



uj°g p ^g p ' k * 
uj°g p ' k g p '" 



1 tu 
1 u 

1 4 w%g p ^g p ^ 



(40) 



Assume that we apply the method to a circular cluster. 
Then we can construct a sub-sample of triplets of galaxies 
inside a ring centered on the cluster center. Each galaxy (i, 



j or k) of the triplet is associated with an angle, 



aj,Oj,afc. 



Hence, we can replace each measured ellipticity (e^, ej 

and 



e,-e 



and 6fe) by the tangential ellipticities (e$ 
£fce 2Mfe ), the above determinant is replaced by: 
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(41) 



We can see that the arguments of the elements of the third 
column are randomly distributed. Therefore the system- 
atic bias vanishes as the inverse of the square root 
of the number of triplets and becomes negligible, whereas 
the main term is not modified. 

Actually the signal to noise ratio of the method (as de- 
fined in (29)) can also be increased in the case of a circu- 
lar geometry because the number of triplets becomes large 



enough to enable a stringent selection within them: cither 
by the angular environment of the three sources (see 5.2), 
either by their redshift (see 5.3.1). 

5.2. Choosing the triplets of galaxies 

This section refers to the necessity of rejecting of the sam- 
ple sources for which another galaxy may play the role of 
a galaxy lens. The angular radius of the circle in which we 
can consider that a galaxy is isolated without important 
local perturbation of another galaxy depends on the real 
mass distribution of galaxy halos. Therefore it is some- 
what difficult to provide accurate estimate of this circle 
with our present-day knowledge of the distant galaxy ha- 
los. As a rough simplification we can consider each galaxy 
as an isothermal sphere and then choose the angular ra- 
dius such that the induced shear (and convergence) is 
lower than 5%. This conservative approximation applies 
not only in the section 4.2.1, but also in 4.2.2, where R gal 
can thus be chosen small (« 2%) in order to decrease the 
systematic bias A^^. 

Of equal importance is the question of the maximum dis- 
tance Ar tolerable between each component of the triplet. 
Or, in the case of a circular cluster analysis, what is the 
thickness of the rings ? To answer we need to balance two 
opposite effects: first, while Ar increases, < > de- 
creases proportionally to about 1 — Ar/r (if the shear of 
the cluster is a law in 1/r, as for an isothermal sphere), 
so the signal to noise ratio decreases; second, in the same 
time the number of triplets increases making possible a ju- 
dicious selection of the redshifts in the triplets (see section 
5.3.1). After realistic simulations (using a number density 
of galaxies equal to 50 galaxies arcmin -2 ), the optimized 
value that has been selected for the first presentation of 
the method is 20 arcseconds separation. In practice, this 
can be optimized accordingly to the shear map of the clus- 
ter. 

5.3. Redshift optimization 

In this section we suppose that the redshift distribution 
n(z) and the variance of the redshift error o-\ z (z) are 
known. We then investigate if the signal to noise ratio 
of the method can be improved thanks to either a choice 
in the redshifts of the triplets or a selection in redshift of 
the clusters. In all the oncoming considerations the num- 
ber density of background galaxies has been taken equal 
to 50 galaxies arcmin -2 ). 

5.3.1. Redshifts in triplets 

From the matrix form (15), it appears obvious that if two 
of the three redshifts of a triplet are very close together 
then the resulting value of G^ k m is nearly zero and thus 
do no bring any cosmological information to G. Therefore 
it is necessary to set two minimum redshift differences 
Ay = zj — Zi and Ajk = Zk — Zj below which the triplets 
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will not be rejected. On the other hand, if these minima 
are too high then too many triplets will be rejected, in- 
creasing the terms 1 / y/N and ac n ■ The right balance is 
obtained for Ay and Ajk greater than about 0.05. 
This small value has been obtained from simulations. It 
proves that in the balance, the second effect is stronger. 
It means that with this density of galaxies, any stringent 
selection makes <jgn decrease rapidly. This will not be 
true anymore if we can select triplets within rings, for 
a nearly circular cluster, or if we use observations with 
NGST (which will increase significantly the number den- 
sity of galaxies). In these two cases , the balance becomes 
favorable to selections and so to an increase of the signal 
to noise ratio of the method. 

The second concern about the redshift selection is the 
difference between Zi and the redshift of the lens, A&. 
Once again, the optimization arises with the balance of 
two competing effects : if Zi is too close to zi then the 
noise is increased by a nearly infinite value of Aw ; ; how- 
ever, in the same time the gradient of G™^ n is greater. 
Our simulations show that the right balance is obtained 
for An greater than about 0.05, and the same discussion 
as above applies. 

5.3.2. Redshifts of clusters 

We can also find the best balance between the last two 
effects when the redshift of the cluster changes. The most 
favorable cluster redshifts are strongly dependent on the 
shape of <ta 2 (z) : the redshift of the cluster must be in 
the area where the errors on the redshifts are as small as 
possible for Auoi to be non dominant in the signal to noise 
ratio. 

Besides, if the redshift of the cluster is too low then 
u>i w uij w ujk whatever the triplet is, the S/N decreases. 
With ctaz(z) approximately equal to 0.03 for redshift infe- 
rior to 0.8 and equal to 0.1 for other redshifts, simulations 
show that clusters at redshift between 0.3 and 0.5 are the 
most favorable to the method. 



6. Simulations and results 

In the simulations, we assumed that the redshift of each 
galaxy was known and we took a redshift distribution n(z) 
which represents both the usual peak of galaxies near 
z = 1 (described by the redshift distribution p(z) from 
Brainerd ct al. 1996) as well as a more distant popula- 
tion of faint blue galaxies suggested by deep spectroscopic 
survey and that seem detectable in the selected color bin 
B = 26 - 28 (see Fort ct al 1996, Broadhurst 98): 



n(z) = aip(z) + a 2 p(z — 2) with 



p(z) 



(3z 2 



(§) ^ 



cxp 



(42) 
(43) 



with z = 1/3 and (3=1, the p(z) and p(z — 2) distribu- 
tions give respectively an average redshift equal to 1 and 
3. With ax = 0.7 and a 2 = 0.3, the fraction of high and 
low redshift galaxies found by Fort et al is respected. It 
concerns a single observation of very faint galaxies in B 
and I color through the cluster C10024+1654. Although 
Broadhurst seems to have confirmed the result of this ob- 
servation with the Keck telescope, the redshift distribution 
may not be exactly representative of the general redshift 
distribution of the faint galaxy population for which we 
can determine good color redshift. However, it is impor- 
tant to remember that the triplet method actually test 
the convexity of the lensing factor curve u{z) (see Fig.l). 
Therefore, the method mostly depends on the total num- 
ber of galaxies above z = 1 that we can detect and for 
which we can determine an accurate color redshift (high 
signal to noise for B photometry is essential). If a deep 
multicolor photometry from U to K gives a broader red- 
shift distribution with more galaxies at intermediate red- 
shift (z = 1.5—2) the result of the triplet method will be 
better since the new distribution will increase the num- 
ber of possible triplets with significant information on the 
cosmology. 

We have also put realistic uncertainties in the photometric 
redshift of galaxies. We assume it is a Gaussian distribu- 
tion with a redshift dispersion equal to 0.03 for redshifts 
lower than 0.8 and as large as 0.1 for redshifts greater than 
0.8 (see Brunner et al. 1997 and Hogg et al. 1998). 
The intrinsic ellipticity distribution is chosen as follows : 



Pe s (es) 



1 



exp 



cr es 



(44) 



with the intrinsic ellipticity dispersion a es — 0.15 (from 
Seitz & Schneider 1997). 

For errors on the measured ellipticities, we chose the same 
distribution, with a dispersion equal to 0.02. 
For the clusters, the redshift have been taken equal to 
0.4, and their projected mass density has the following 
analytical shape : 



K (r,0) = 



rVl + e cos 26*' 



(45) 



The cluster ellipticities were chosen randomly between 
and 0.5. It is unimportant to include a core radius since 
its influence in the weak lensing area is negligible (Fur- 
thermore, we recall that the principle of the method is to 
cancel any dependency on the mass map of the cluster). 
k is taken to correspond to a sample of high velocity dis- 
persion (1000 km/s). 

The size of the observation window has been taken simi- 
lar to what will give the VLT instrument FORS : 6' x 6'. 
Within this window, a circle of 90" radius representing the 
arclets and strong lensing regimes was not considered. 
So, with an average of about 50 galaxies arcmin" 2 , each 



12 L.Gautret 1 , B. Fort 2 , Y.Mellier 1,2 : Constraining (SI, A) from weak lensing in clusters: the triplet method. 



VLT observation of a cluster contains about 1000 back- 
ground galaxies. 

Now, the question that the simulations have to solve is: 
what is the accuracy we can reach on the cosmological pa- 
rameters for a given number of observed clusters ? 
Figure (4) gives the contours pn^f^A) for simulations 
concerning 1000 clusters in three cases : first case, Sl = 
0.3 and A Q = 0.7; second case, Sl = 0.3 and A c = 0; third 
case, Sl a = 1 and A Q = . Figure (5) gives the same for 
100 clusters. 

The results of these simulations are promising. They prove 
that with about 100 clusters A^q^ (in the case Sl + A = 1) 
and fl±oi° 5 (in the case A = 0) can be reached (with a 
70% confidence level). 100 clusters correspond to about 
20 nights of VLT observation. SI = 0.3 and St = 1 uni- 
verses can also be separated at a 2a = 95% confidence 
level with the same time of observation: about 20 VLT 
nights. Therefore, even a modest observing campaign on 
a VLT could provide interesting constraints on A. 



7. Conclusion 

We have presented a method to constrain the cosmologi- 
cal parameters (SI, A). It uses gravitational distortion pro- 
duced by weak lensing around clusters and the photo- 
metric rcdshift of lensed galaxies. This purely geometri- 
cal method is insensitive to the lens modeling and can 
be directly applied to real data, that is the ellipticities of 
galaxies as observed from optical images. We have calcu- 
lated the statistical noise coming mainly from the intrinsic 
source ellipticity and performed realistic simulations. 
The main result is that with a short program of obser- 
vations (about 20 VLT nights) a constraint could be pro- 
vided on the value of the cosmological constant. Using the 
observation of 100 clusters, we can reach A+q^ m the case 
O + A = 1 or SIIq^ in the case A = (at a la level). 
These results could be even refined down to 10% accuracy 
on (SI, A) with the use of NGST. Indeed, the NGST looks 
perfectly suited for this method since it increases the num- 
ber density of observed background galaxies. 
We have estimated the amplitude of the systematics due to 
the presence of background structures with a multi-lcnsing 
model. It turns out that the shift amplitudes on the (SI, A) 
determination are about 0.05. One part of this systematic 
can be directly corrected. It concerns the perturbative ef- 
fect due to galaxy-galaxy lensing (the correction of this ef- 
fect need an approximate mean potential of galaxies and 
will be provided by the incoming galaxy-galaxy lensing 
investigations) . It also concerns the effect of background 
structures integrated along the line of sight (the correction 
of this effect requires calculations that takes in account the 
non linear evolution of large scale structures). It could be 
validated by ray tracing simulations. 
In conclusion, we have shown that the systematic effects 
could be very well controlled by a judicious selection cri- 
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Fig. 4. Contours of the (SI, A) probability distribution ob- 
tained from a simulation on 1000 clusters in the cases 
(SI, A) = (0.3, 0.7) (top panel) , (SI, A) = (0.3, 0) (middle panel) 
and (1,0) (bottom panel). We give the la — 68% (grey) and 
2a — 95% (dark) confidence levels. 



L.Gautret 1 , B. Fort 2 , Y.Mellier 1,2 : Constraining (f2, A) from weak lensing in clusters: the triplet method. 13 



P„(Q,A) = 0.3 A=0.7 




P„(Q,A) Q = 0.3 A=0 




0.0 0.2 0.4 0.6 0.8 1.0 




Fig. 5. Contours of the (fi,A) probability distribution ob- 
tained from a simulation on 100 clusters in the cases 
(fi, A) = (0.3,0.7) (top panel) , (fi,A) = (0.3,0) (middle panel) 
and (1,0) (bottom panel). We give the la — 68% (grey) and 
2a = 95% (dark) confidence levels. 



teria of the clusters and the lensed galaxies of each triplet. 
The (f2, A) degeneracy of the triplet method is orthogonal 
to the one of the classical (m, z) of supernovae searches. In 
this regards, when combined to the supernovae approach, 
the triplet method using VLT or NGST data could be 
extremely efficient. Therefore it seems important to inves- 
tigate more deeply the possibility to use cluster lenses as 
practical tests to constrain the curvature of the Universe. 
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Appendix A : The observed ellipticity of a back- 
ground source 

This appendix has two goals: firstly to recall the practi- 
cal way to calculate the image ellipticity ej — eie 2te ' of 
a source from the image of a galaxy; secondly to express 
this ellipticity as a function of the potential of the lens- 
ing cluster through the local convergence n and complex 
shear 7 = -fe 2l0L (as it will be showed, the dependency in 
k, 7 holds in the complex term g = jz^) and the intrinsic 
source ellipticity eg = ese 2l0s . 

For each background galaxy we can calculate a second or- 
der momentum matrix M 1 , either directly from the image 
of the galaxy (see Bonnet 1995) either from the ACF of 
the single galaxy (see Van Waerbeke et al. 1997). Then the 
complex image ellipticity is derived from the momentum 
matrix through (see Seitz & Schneider 1997): 



e/ 



M li - M 22 + 2iM{ 2 
M li + M k + 2 y /det(M I ) ' 



(46) 



It is easy to check that this definition stays in agree- 
ment with equation (5). We also use a second order mo- 
mentum matrix M s for the source as it would be seen 
with no lens. 

The effect of the lens on the ellipticity of the background 
galaxy can then be summarized in the matrix equation 
(see Bonnet 1995): 



M 1 = 
A = 



AM S l A 



with 



\A\ 2 



(1-k) 2 - 7 2 ' 
In what follows and above, we use the notations 



T 29 



<s(26) -sin(20) \ 



v sin(26») cos(26>) J ' 

2g /cos(20) sin(2fl) \ 
^sin(26>) -cos(20) J ' 

Equation (46) leads to 

(1° + gJ 20 )I Cp J ) + gJ 20 ) 

rl _ \ U 1+e / 



(47) 
(48) 

(49) 
(50) 



M 1 = 



(51) 



|l-5 2 l 

Using the following properties : 

jaj fj = jaj = ja -f3^ jaj/3 = ja +f3 _ ^ 

equation (51) gives: 

|l- fl 2 |(l-e|)M^= (l + e 2 )(l+g 2 ) 

+4ge s cos(2(9 s - 6 L ))]I° 
+2e s J 20s +2g(l + e 2 )J 2e - [b6) 
+2g 2 e s J i0L - 20s . 



Finally, whatever the value of g is (meaningly either in 
the weak or strong lensing regime) the measured image 
ellipticity writes: 



ei 



e s + (1 + e|)g + e^g 2 



max(l, g 2 ) + e 2 s min{l, g 2 ) + e* s g + e$g* ' 



(54) 



which, in the weak lensing and arclet regimes (g < 1) 
simplifies into (to the third order in es and g): 



£/ = (1 -.9 2 )es +9 - 



* 2 
9 e s , 



(55) 



where the term g*e- 2 s is negligible regarding to the term 



es and vanishes in 1/vJV (because the argument of this 
term behaves randomly) in the operator G as the different 
noises described in section 3.4. Therefore the final equa- 
tion for the image ellipticity writes : 



e/ = (1 - g 2 )e s +g. 



(56) 



It is worht noticing here that the term (1 — g 2 ) is 
in favor of the method developed in this paper. Indeed, 
with the conditions used in the simulations (section 6), 
the mean value of this term is about 0.85. It means that 
the noise coming from the intrinsic source ellipticity is 
lowered by 15% or, for a given signal to noise ratio, the 
number of required clusters is lowered by 28%. 

Appendix B : The case of a perturbing lens 

This section will study the influence of a perturbative lens 
(galaxy, cluster or higher scale structure) located behind 
the principal lensing cluster on the measured ellipticity of 
a background source. 

The calculation of the total amplification matrix in the 
case of two (or more) consecutive lenses has already been 
done (see for example Kovner 1997). The result can be 
simply noticed as follows: 



A- 1 = 
L = 
L p = 



kI° + 7 J 



K p I° +1 p J 20 



L + cL^L with 

20 - and 

p 



(57) 
(58) 
(59) 



The p upper index refers to the perturbative lens; k p 
is defined as k in equation (3), just changing the principal 
lens noticed I by the perturbative lens l p for the definition 
of the perturbative critical surface mass density T, p rit . The 
c coupling factor is c = ®ll pD °s 



Equation (57) rewrites: 



A' 1 = (1 



+ CKK P )I° - 7(1 - CK P )J 



20, 



(60) 



- 7 P (1 - cn)J 20 ? + c 7 7 P / 2 (^-H 
which can be inverted into 
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A = 



[(1 — K — K T CKK }I + 7 (1 - CK F ) J Wl (61) 
+ 7 P (1 - CK)J 20 r + m P I J ( ,l - ,f )]. 



Contrary to the case of a single lens where the Am- 
plification matrix is symmetrical, here appears an anti- 
symmetric term, the last one of the above equation. How- 
ever, even if the Amplification matrix is non symmet- 
ric, it keeps the propriety to transform an ellipse into 
another ellipse from the source plan to the image plan. 
Indeed the transformation from the source to the image 
can be represented by the vector notation X$ = A~ x Xj. 
To say that the source is an ellipse is equivalent to say 
that exists a matrix E such that t X t s EEXs — 1 (it comes 
from the fact that every positive symmetric matrix can be 
written as t EE). So the position of image vectors writes 
t X\EEXi = 1 with E = EA^ 1 , which proves that the 
image is another ellipse. We thus can search for the im- 
age ellipticity of a source distorted by both lenses, as a 
function of the source ellipticity and the potential of the 
principal and perturbative lenses. We proceed the same 
way as in annex A. 

In the following, we only keep second order terms in 
k, k p , 7, 7 P and es- So we can use the definitions: 



2. Secondly a scalar correction (1 — c)k p on the cosmo- 
logical term g, due to coupling between the main and 
the perturbative lenses. 

The simple form of equation (66) can be easily gener- 
alized to the case of a succession of many lenses noticed 
L° (instead of the principal lens), L 1 (instead of the per- 
turbative lens), L 2 , L n between the observer and the 
source noticed S. If the coupling factor between the lenses 
i and j (with Zi < z J ) is noticed cfj = jy^-jff^ , then 

equation (55) giving the amplification matrix can be gen- 
eralized into: 



)(J° - V) = (1° - L n ) 



(I°-L°), (67) 



where the ® symbol defines a particular matrix prod- 
uct with the following properties: 7° ® L l — L l ® 1° = L l , 

U ® L l — " S T 3 Ti r.fc (Qi 13 i<* Ti — *S „S rk j 



I - J I' etc... 



We can remark that in the case of two very close lenses, 
their associated coupling factor c is zero. With the same 
calculations and approximations as above we obtain: 



g p 



= (1-k-(1-c)k p Y 
= (1 - (1 -c)n-K P y 



7 



and rewrite the amplification matrix: 

1° + gJ 20 t + g p J 20 r + c~gg P I 2{eL - 6p) 



A = 



9 2 



xP2 



2(1 - c)Re(gg p *) 



(62) 
(63) 



(64) 



Finally, with a calculation similar to the one done in annex 
A, we obtain for the image ellipticity: 



e/ = 



(l-\g + g p \ 2 )e s + g + g p 



-4cRe(gg p *(g + g p ) 
which we simplify into 



2cg*g p e s 



e I = (l-\g + g p \ 2 )e s +g + g 



(65) 

4(9* + 9 p n 

(66) 



ej = (l-|^| 2 )e s - 

j=0 

A* - (1-^-^(1- 



i—n 
i=0 



with 



(68) 
(69) 



In this last equation we have cancelled all the terms 
higher to the second order and oriented randomly (we con- 
sider that the source ellipticity, the principal and pertur- 
bative shears are independent. Terms higher to the second 
order with a random orientation give a negligible noise - 
regarding to the noise coming from the intrinsic source 
ellipticity- in the calculation of G). 

Comparing equations (62) and (66) we can see the effect 
of the perturbative lens on the measured ellipticity: 

1. Firstly a complex additive correction g p which orien- 
tation is non correlated to the one of the cosmological 
term. 
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